Hubbard model on an anisotropic checkerboard lattice at finite temperatures: Magnetic and metal–insulator transitions
Liu Hai-Di
State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China

 

† Corresponding author. E-mail: hdliu@wipm.ac.cn

Abstract

We study magnetic and Mott transitions of the Hubbard model on the geometrically frustrated anisotropic checkerboard lattice at half filling using cellular dynamical mean-field theory. Phase diagrams over a wide area of the parameter space are obtained by varying the interparticle interaction strength, geometric frustration strength, and temperature. Our results show that frustration and thermal fluctuations play a competing role against the interactions and in general favor a metallic phase without antiferromagnetic order. Due to their interplay, the system exhibits competition between antiferromagnetic insulator, antiferromagnetic metal, paramagnetic insulator, and paramagnetic metal phases in the intermediate-interaction regime. In the strong-interaction limit, which reduces to the Heisenberg model, our result is consistent with previous studies.

1. Introduction

Geometrical frustration in the strongly correlated systems has attracted considerable attention recently and was widely studied.[19] Strong interparticle interactions are known to play a key role in many interesting phenomena in condensed matter physics, including Mott transition and high-temperature superconductivity,[1012] but it still remains a challenge to fully understand them due to the break down of perturbative approaches. Geometric frustration, on the other hand, describes a phenomenon where, due to geometry, global minimization of inter-particle interaction energies cannot be trivially achieved, resulting in various non-trivial ground state configurations, such as quantum spin liquid, spin ice, and spin glass.[1315] Antiferromagnetically interacting spins arranged with certain geometries, e.g., on a triangular lattice, are a typical type of systems that host geometric frustration. Geometric frustration also plays an important role in interacting electronic systems, especially in half-filled Hubbard models on geometrically frustrated lattices, whose strongly-interacting limit reduces to spin models. It is intriguing to investigate the competition between interaction and geometrical frustration in such systems, which exhibit fascinating interplay between magnetic transitions and metal–insulator transitions, leading to novel phases and rich phase diagrams.

For antiferromagnetically interacting spins, the simplest two-dimensional (2D) geometry that leads to geometric frustration is a triangle, which can be extended to a triangular lattice or a kagome lattice. In three-dimension (3D), the simplest example is a tetrahedron, whose 3D extension describes magnetic properties of some transition-metal oxides, such as LiV2O4 and Tl2Ru2O7.[57] Similarly, one can construct a 2D lattice with a layer of vertex-sharing tetrahedra. It is topologically equivalent to a 2D square lattice with selected diagonal links, as shown in Fig. 1(a), which is called a checkerboard lattice.[2,1630] There are two sites in its unit cell, corresponding to its two sublattices, which we label with A and B, and its first Brillouin zone is show in Fig. 1(b). For a Hubbard model, the linking strength is characterized by the hopping energy along the link. In general, the non-diagonal and diagonal links can have different strengths, which correspond to hopping energies t1 and t2, respectively. When , it reduces to a 2D square lattice which is frustration free. Therefore, one can use t2 as a measure of frustration and study the effect of geometric frustration by varying it. The special case where is called the fully-frustrated case, or isotropic in the sense that it resembles a 2D layer of vertex-sharing regular tetrahedra. There have been lots of works on antiferromagnetic Heisenberg models on checkerboard lattices.[1926] The ground state is shown to have macroscopic degeneracy without magnetic orders.[19] Previous studies on Hubbard models on checkerboard lattices mainly focused on the isotropic case. A study based on renormalization group method and mean-field analysis shows a (semi)metal–insulator transition at half filling.[2] Reference [28] presents a path-integral renormalization group analysis, and finds a first-order phase transition to a plaquette-singlet insulating phase at zero temperature. In Ref. [29], based on an auxiliary field decomposition, the authors applied Monte Carlo and variational methods, and obtained a low-temperature phase diagram, which contains exotic phases like an insulating phase with flux-like correlations and a 120° correlated state. The path-integral renormalization group method has also been applied to the anisotropic case where to obtain a zero-temperature phase diagram of the system.[30] Compared with the isotropic case, the anisotropic one is less understood and demands further study. In this article, we take the half-filled Hubbard model on the anisotropic checkerboard lattice at finite temperatures to study the interplay between geometrical frustration, interparticle interaction, and thermal fluctuations.

Fig. 1. (a) The checkerboard lattice. The blue (red) dots represent lattice sites of the sublattice. Black solid and black dashed lines represent nearest-neighbor and diagonal links, respectively. We choose the distance between the nearest neighbors as the length unit in this article. The yellow box containing four lattice sites marks one cluster used in the CDMFT calculation. (b) The first Brillouin zone (BZ) of the checkerboard lattice. We also show the reduced first Brillouin zone for the chosen cluster in CDMFT calculations. High-symmetry points of the BZ are marked with , K, and M. (c) Non-interacting band structure of a tight-binding model on a homogeneous ( ) checkerboard lattice. (d) Non-interacting density of states when .
2. Model and method

We consider the standard Hubbard model on the anisotropic checkerboard lattice at half filling. Its Hamiltonian can be written as

where ( ) is the creation (annihilation) operator of fermions on site i with spin σ, is the number operator, t1 (t2) is the nearest-neighbor (diagonal) hopping energy, and U is the strength of the on-site interaction. ( ) runs over the nearest-neighbor (diagonal) links on the lattice as shown in Fig. 1(a). We choose the distance between the nearest sites as the length unit throughout the article.

To start with, we first look at the non-interacting (U = 0) band structure, which is readily obtained by diagonalizing the kinetic-energy terms (the first two terms on the right hand side of Eq. (1)). The energy of the resulting two bands is

We plot the dispersions along high symmetry lines of the first Brillouin zone (see Fig. 1(b)) with different frustration strengths (t2) at half filling in Fig. 2. When , there is no frustration and the lattice effectively reduces to a square lattice with perfect Fermi-surface nesting, which leads to a Mott insulating phase with any finite positive U at zero temperature. For the fully-frustrated case or isotropic case with , the upper band becomes completely flat; the lower dispersive band touches it quadratically at the corner of the first Brillouin zone. We also plot the band structure over the whole first Brillouin zone and the corresponding density of states (DOS) of the homogeneous case ( ) in Figs. 1(c) and 1(d), respectively. At half-filling, the chemical potential sits at this band-touching point or right at the flat band, which causes divergences in perturbative calculations.[2,3032] In this article, we take a numerical approach to this model, covering a wide range of frustration strength from the frustration-free case to the fully-frustrated one at finite temperatures.

Fig. 2. Non-interacting band dispersions along high symmetry lines (see Fig. 1(b) with different frustration strengths at half filling. We choose as the energy unit. is the chemical potential at half filling.

Our analysis of the interplay between geometrical frustration and electron–electron interaction in the half-filled checkerboard lattice is bases on the cellular dynamical mean-field theory (CDMFT)[3335] calculation combined with continuous time quantum Monte Carlo method (CTQMC),[36,37] which has been widely applied to strongly correlated geometrically frustrated systems. The main idea of CDMFT is to map a system on a lattice to an impurity model by considering a cluster embedded in a self-consistent effective Weiss field. This impurity model is then solved using CTQMC, which updates the effective Weiss field and completes one self-consistent loop. Using this CDMFT approach, we calculate the single-particle Greenʼs function at imaginary frequencies, which is then analytically continued to real frequencies via the maximum entropy method.[38] With the single-particle Greenʼs function, one can further calculate physical quantities like DOS, magnetization, and spectral functions.

In our calculation, we choose four sites marked by a yellow circle as one cluster as shown in Fig. 1(a), and the corresponding reduced first Brillouin zone is shown in Fig. 1(b). We set as the energy unit so that t2 can be regarded as the frustration strength or the anisotropic parameter. Convergence is achieved by iterating the CDMFT self-consistent loop, and in each iteration 107 QMC steps are performed. We work with the grand canonical ensemble and the half-filling condition is realized by adjusting the chemical potential.

3. CDMFT results

We focus on two important aspects of phase transitions in this system, metal–insulator transitions and magnetic transitions. The former is identified by looking at the DOS, where an energy gap at the chemical potential indicates an insulating phase, and vice versa. The latter is determined through the corresponding magnetic order parameter.

3.1. Metal–insulator transitions

For fixed frustration strength t2, the system goes through the typical Mott transition to a insulating phase when the interaction strength U is increased. This can be seen from Fig. 3, where each subfigure corresponds to a fixed frustration and shows plots of DOS at different interaction strengths.

Fig. 3. Density of states at different interaction strengths with (a) T=0.05, , (b) T=0.05, , (c) T=0.05, , and (d) T=0.1, . The frequency ω (which is also the energy, because we set the reduced Planck constant to unity) is counted from the chemical potential.

For example, Figure 3(a) is the non-frustrated case with temperature and . One can see that when U is small, there is a Fermi-liquid-like peak near the Fermi energy, which is then split into a pseudogap structure when the interaction increases.[39,40] When the interaction strength is larger than a critical value, which is for this case, the DOS develops an obvious energy gap at the chemical potential, which indicates that the system undergoes a Mott transition from a metallic phase to a Mott insulating phase. Similarly, from Figs. 3(b) and 3(c), one can read out and 4.3 for and 0.5, respectively. We can see that the critical interaction strength of Mott transition increases with the frustration strength increasing. We also calculated the DOS with a higher temperature and which is plotted in Fig. 3(d) and shows . It is consistent with the general picture that thermal fluctuations make it more difficult to drive a system to an ordered phase.

To gain a finer view of the Mott transition, it is instructive to look at the double occupancy (Docc) defined as , where is the number of sites in one cluster. It indicates the average possibility of two particles occupying the same site. We calculated Docc as a function of U at different temperatures and frustration strengths, and the result is shown in Fig. 4, where the arrows with different colors mark the corresponding phase transition points. It can be seen that Docc decreases with an increasing interaction strength, which is a signature of Mott transition. In the large U limit in the Mott insulating phase, the system is almost singly occupied, i.e., most particles are localized and hopping between lattice sites is forbidden by the strong on-site interactions. Moreover, when temperature is high (see plots with T=0.1 and 0.2 in Fig. 4), Docc decreases smoothly as U increases, indicating that the Mott transition is continuous. At lower temperature (see plots with T=0.05 in Fig. 4), however, there is a visible discontinuity in Docc around Uc. This suggests that the transition is a first-order transition at low temperatures.[9,28,41]

Fig. 4. Double occupancy as a function of the interaction strength U at different temperatures and frustration strengths. Arrows with different colors mark the corresponding positions of the critical U of the Mott transition.

In addition, we further calculate the distribution of the spectral weight at the chemical potential with different interaction strength. It is a good characterization of interaction effects and can be directly accessed via angle resolved photoemission spectroscopy (ARPES) experiments. The result with T=0.05, is shown in Fig. 5. When U = 3, smaller than the critical interaction strength for Mott transition, the spectral function has sharp peaks, indicating well-defined quasiparticles. Since it is calculated at zero frequency, it also shows a clear Fermi surface structure.[8] In the Mott insulating phase when U = 6, sharp peaks disappear and no clear Fermi surface is shown.

Fig. 5. Distribution of the spectral weight at zero energy (the chemical potential) with (a) U = 3 and (b) U = 6. We choose T=0.05 and . For easy comparison, we use the same scale for spectral weight in these two plots.
3.2. Magnetic transitions

For magnetic transition, we focus on the antiferromagnetic order and define the order parameter as , where if i belongs to sublattice (see Fig. 1(a)). The calculation of m is straightforward with the single-particle Greenʼs function obtained. We plot m as a function of U at a fixed temperature and frustration strength (T=0.05, ) together with the energy gap obtained from DOS under the same condition in Fig. 6. When the interaction is weak, the magnetic order parameter is zero and the system is in a paramagnetic phase. At U = 4.2, the system develops nonzero m with opposite signs between A and B sublattices, indicating an antiferromagnetic order. This is consistent with the fact that the strong interacting limit can be described by an antiferromagnetic Heisenberg model with frustration induced by t2.

Fig. 6. Antiferromagnetic order parameter m and single-particle energy-gap as a function of the interaction strength U when T=0.05 and .

With the information of the energy gap, one can further see that the development of a nonzero energy gap and a nonzero antiferromagnetic order do not occur simultaneously when T=0.05, . This leads to three phases, the paramagnetic metal phase (PM) when , the antiferromagnetic metal (AFM) phase when , and the antiferromagnetic insulator phase (AFI) when .

3.3. Phase diagrams

Based on the above analysis of Mott transitions and magnetic transitions, we are ready to present phase diagrams which cover a wide area of the parameter space.

We first look at the result near zero temperature. Because zero temperature is not accessible through CTQMC, we approach the zero-temperature case by calculating at low temperatures, and a low-temperature (T=0.05) phase diagram in the t2U plane is shown in Fig. 7. It clearly shows the competing effect of magnetic frustrations and inter-particle interactions. An increasing interaction strength drives the system into an insulating phase with an antiferromagnetic order, while frustrations favor a metallic phase and suppress the magnetic order. Furthermore, the critical t2 where antiferromagnetic order disappears tends to saturate when U increases. This suggests that when the frustration exceeds a critical value (here at ), it completely destroys the anti-ferromagnetic order even at the large U limit. This observation is consistent with previous studies that focus on the effect of frustration and deal with the Heisenberg model,[26] which effectively describes the large U limit of the model we considered here. The critical value of also agrees well with previous studies using different approaches.[30] In addition, there also exists a narrow window for an antiferromagnetic metal phase between the paramagnetic metal and antiferromagnetic insulator phases.

Fig. 7. Phase diagram in the t2U space with T=0.05. The red line is the phase boundary between the paramagnetic phase and the antiferromagnetic phase. The black line indicates the Mott transition (dashed line for first-order transitions and solid line for continuous ones). These two kinds of transitions divide the phase diagram to regions corresponding to paramagnetic metal (PM), antiferromagnetic metal (AFM), antiferromagnetic insulator (AFI), and paramagnetic insulator (PI) phases. The red star marks the critical point of the magnetic transition of the Heisenberg limit ( ) obtained by previous studies.[26] It can be seen that our calculation agrees with it well.

We also investigated the effect of thermal fluctuations by presenting a phase diagram on the TU plane at a fixed frustration strength , as illustrated in Fig. 8. One can see that thermal fluctuations and interparticle interactions also play competing roles in determining the phases of the system. As expected, when the temperature increases, the thermal fluctuations eventually destroy antiferromagnetic order as well as the Mott insulating phase and drive the system to a paramagnetic metal phase. During this transition, there is also an antiferromagnetic metallic phase in between.

Fig. 8. Phase diagram in the TU space at .
4. Discussion and summary

To summarize, we have studied Mott transitions and antiferromagnetic transitions on a half-filled anisotropic checkerboard lattice by means of the cellular dynamical mean field theory combined with the continuous-time quantum Monte Carlo method. We benchmark our calculation by looking at the critical frustration strength in the large-U limit, which can be effectively described by an antiferromagnetic Heisenberg model, and our result is consistent with previous studies. The interplay between frustrations (t2) (or anisotropy with the isotropic case being most frustrated), interparticle interactions (U), and thermal fluctuations (T) are investigated by varying the corresponding parameters. Through numerical simulations in a wide area of the parameter space, we extract phase diagrams of the system in the t2U and TU spaces. At fixed temperature and frustration strength, the interaction drives a Mott transition. When the interaction is sufficiently strong, the lattice sites are nearly singly occupied, and the system is in the insulating phase. This transition is accompanied by a magnetic transition, which tends to produce antiferromagnetic order when the interaction is strong. Thermal fluctuations and frustrations, on the other hand, play a competing role and favor metallic phases, and the magnetic order can be completely suppressed when t2 and T exceed certain values. Besides, they further complicate the phase diagram by opening a narrow window for an antiferromagnetic metal phase. We hope that our study can provide useful information for future research on strongly-interacting systems on checkerboard lattices and related materials.

Reference
[1] Hagemann I S Huang Q Gao X P A Ramirez A P Cava R J 2001 Phys. Rev. Lett. 86 894
[2] Fujimoto S 2002 Phys. Rev. Lett. 89 226402
[3] Kimura T Kuroki K Arita R Aoki H 2004 Phys. Rev. B 69 054501
[4] Ohashi T Momoi T Tsunetsugu H Kawakami N 2008 Phys. Rev. Lett. 100 076402
[5] Takeda T Nagata M Kobayashi H Kanno R Kawamoto Y Takano M Kamiyama T Izumi F Sleight A W 1998 J. Solid State Chem. 140 182
[6] Sakai H Kato M Yoshimura K Kosuge K 2002 J. Phys. Soc. Jpn. 71 422
[7] Kondo S Johnston D C Swenson C A Borsa F Mahajan A V Miller L L Gu T Goldman A I Maple M B Gajewski D A Freeman E J Dilley N R Dickey R P Merrin J Kojima K Luke G M Uemura Y J Chmaissem O Jorgensen J D 1997 Phys. Rev. Lett. 78 3729
[8] Parcollet O Biroli G Kotliar G 2004 Phys. Rev. Lett. 92 226402
[9] Ohashi T Kawakami N Tsunetsugu H 2006 Phys. Rev. Lett. 97 066401
[10] Imada M Fujimori A Tokura Y 1998 Rev. Mod. Phys. 70 1039
[11] Dagotto E 1994 Rev. Mod. Phys. 66 763
[12] Rohringer G Hafermann H Toschi A Katanin A A Antipov A E Katsnelson M I Lichtenstein A I Rubtsov A N Held K 2018 Rev. Mod. Phys. 90 025003
[13] Shimizu Y Miyagawa K Kanoda K Maesato M Saito G 2003 Phys. Rev. Lett. 91 107001
[14] Binder K Young A P 1986 Rev. Mod. Phys. 58 801
[15] Nisoli C Moessner R Schiffer P 2013 Rev. Mod. Phys. 85 1473
[16] Xu Y Xiong Z Wu H Q Yao D X 2019 Phys. Rev. B 99 085112
[17] Zeng T S Zhu W Sheng D 2018 npj Quantum Materials 3 49
[18] Wu H Q He Y Y Fang C Meng Z Y Lu Z Y 2016 Phys. Rev. Lett. 117 066403
[19] Moessner R Chalker J T 1998 Phys. Rev. B 58 12049
[20] Berg E Altman E Auerbach A 2003 Phys. Rev. Lett. 90 147204
[21] Moukouri S 2008 Phys. Rev. B 77 052408
[22] Lieb E H Schupp P 1999 Phys. Rev. Lett. 83 5362
[23] Berg E Altman E Auerbach A 2004 Phys. Rev. Lett. 90 147204
[24] Capponi B 2017 Phys. Rev. B 95 014420
[25] Morita K Shibata N 2016 Phys. Rev. B 94 140404 R
[26] Brenig W Grzeschik M 2004 Phys. Rev. B 69 064420
[27] Pollmann F Betouras J J Shtengel K Fulde P 2006 Phys. Rev. Lett. 97 170407
[28] Yoshioka T Koga A Kawakami N 2008 J. Phys. Soc. Jpn. 77 104702
[29] Swain N Majumdar P 2017 J. Phys. Condens. Matter. 29 085603
[30] Yoshioka T Koga A Kawakami N 2008 Phys. Rev. B 78 165113
[31] Isoda M Mori S 2000 J. Phys. Soc. Jpn. 69 1509
[32] Fujimoto S 2001 Phys. Rev. B 64 085102
[33] Kotliar G Savrasov S Y Palsson G Biroli G 2001 Phys. Rev. Lett. 87 186401
[34] Maier T Jarrell M Pruschke T Hettler M H 2005 Rev. Mod. Phys. 77 1027
[35] Park H Haule K Kotliar G 2008 Phys. Rev. Lett. 101 186403
[36] Gull E Millis A J Lichtenstein A I Rubtsov A N Troyer M Werner P 2011 Rev. Mod. Phys. 83 349
[37] Rubtsov A N Savkin V V Lichtenstein A I 2005 Phys. Rev. B 72 035122
[38] Jarrell M Gubernatis J E 1996 Phys. Rep. 269 133
[39] Imai Y Kawakami N 2002 Phys. Rev. B 65 233103
[40] Huscroft C Jarrell M Maier T Moukouri S Tahvildarzadeh A N 2001 Phys. Rev. Lett. 86 139
[41] Yoshioka T Koga A Kawakami N 2009 Phys. Rev. Lett. 103 036401